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MODELLING  SHOCK  INITIATION  AND  DETONATION 
IN  THE  NON  IDEAL  EXPLOSIVE  PBXW-1 15 

David  L.  Kennedy 
ICI  Australia  Operations  Pty.  Ltd. 

PO  Bon  196.  Kurri  Kurri  NSW  2327.  Australia 

David  A.  Jones 

DSTO  Materials  Research  Laboratory 
PO  Box  50,  Ascot  Vale  3032.  Australia 

We  analyse  the  detonics  of  the  non-ideal  explosive  PBXW-115  (also  called  PBXN-111). 
Two  chemical  equilibrium  codes  were  used  to  predict  its  ideal  CJ  state,  with  estimates  of 
its  ideal  detonation  velocity  differing  by  over  1.3  mm.ps1.  A  small  divergent  detonation 
theory  was  calibrated  to  unconfined  detonation  velocity  measurements,  and  used  to 
describe  the  CJ  state  at  different  charge  diameters.  It  was  predicted  that  the  detonation  was 
being  supported  by  about  15%  reaction  at  the  critical  diameter,  with  the  velocity  about 
1.7  mm.ps'1  below  the  ideal  value.  A  finite  element  hydrocode  was  used  to  simulate  a 
variety  of  initiation  and  detonation  tests,  with  the  results  generally  in  excellent  agreement 
with  the  experimental  data.  It  was  found  that  many  experimental  techniques  (including 
those  for  the  measurement  of  ideal  detonation  velocity.  CJ  zone  length,  and  Pop-plots) 
require  substantially  different  interpretation  for  non-ideal  explosives.  _ 


INTRODUCTION 

Composite  explosives  are  used  extensively  in  both 
commercial  and  military  applications.  The  partial  or  complete 
physical  separation  of  their  oxidiser  and  fuel  phases  results  in 
mass  diffusion  reducing  their  reaction  rates  and  increasing 
their  reaction  zone  lengths  in  comparison  with  mono- 
molecular  explosives.  As  a  consequence,  they  exhibit  non¬ 
ideal  performance,  where  their  detonation  velocities  can  be 
considerably  lower  than  the  ideal  value  predicted  by  equilib¬ 
rium  thermodynamic  calculations.  They  typically  have  large 
critical  diameters.  These  factors  all  contribute  to  increased 
difficulty  and  expense  when  performing  characterisation 
experiments,  as  the  charge  weights  involved  are  necessarily 
larger  than  those  involving  ideal  explosives.  This  raises  the 
question  as  to  how  relevant  existing  experimental  techniques 
are  for  the  investigation  of  the  detonics  of  such  explosives. 

This  paper  will  present  a  detailed  evaluation  of  the 
experimental  techniques  applied  to  the  non-ideal  explosive 
PBXW-115  (also  known  as  PBXN-111).  These  include 
detonation  velocity  and  critical  diameter  for  both  unconfined 
and  confined  charges,  comer  turning  ability,  and  shock  front 
curvature,  all  obtained  at  NSWC  by  Forbes  and  co workers1-1 
for  die  US  variant.  Bocksteiner  et  al.3  at  MRU  have  measured 
detonation  velocity  and  critical  diameter  for  both  unconfined 
and  confined  charges  of  the  Australian  variant,  examining  the 
effects  of  RDX  and  A1  panicle  size.  Held4  has  apptied  a 
small-sample  technique  to  the  measurement  of  the  ideal 
detonation  velocity  of  the  German  variant. 

The  results  will  include  a  description  of  the  ideal 
Chapman-Jougct  (CJ)  detonation  state  as  provided  by  two 
equilibrium  thermodynamic  codes,  an  analysis  of  steady 
detonation  in  axisymmetric  geometry,  and  time-resolved 
numerical  simulations  using  the  hydrocode  DYNA2D.5 


IDEAL  DETONATION  —  THE  IDEAL  CJ  STATE 

The  nominal  density  of  PBXW-1 15  is  about  1.79  g.cm  3, 
corresponding  to  voidage  levels  of  less  than  2%.  It  is  cast 
cured,  with  its  nominal  composition  being  20%  cydo- 
irimethylene  trinitrate  (RDX).  43%  ammonium  perchlorate 
(AP),  25%  aluminium  (Al),  and  12%  hydroxy  terminated 
polybutadiene  (HTPB)  as  binder. 

Two  chemical  equilibrium  codes  were  used  to  calculate 
its  CJ  state  for  an  ideal  detonation,  with  the  results  presented 
in  Table  1.  The  first  was  BKW,6  based  on  the  Becker- 
Kistiakowsky-Wilson  equation  of  state  (EoS)  for  the  gaseous 
products  and  the  Cowan  EoS  for  the  solid  products.  The 
second  was  the  ICI  code  IDeX  (standing  for  Ideal  detonation 
of  explosives),  using  an  intermolecular  EoS’  for  the  gaseous 
products  and  the  Mumaghan  EoS  for  the  solid  products. 

TABLE  1 .  IDEAL  CJ  STATE  OF  PBXW- 1 15 


Parameter 

BKW 

IDeX 

Detonation  velocity  (mm.ps1) 

8.010 

6.665 

Detonation  pressure  (GPa) 

27.51 

22.53 

Detonation  temperature  (K) 

5175 

5295 

Heat  of  reaction  (MJ.kg"1) 

6.328 

8.396 

CJ  Gamma 

3.175 

2.534 

The  predictions  from  the  two  codes  differ  to  a  greater 
extent  than  has  been  observed  for  any  other  explosive.  Some 
of  this  would  be  due  to  the  treatment  of  the  binder  —  the 
BKW  calculation  used  an  existing  plasticiser  (C41H74O6) 
from  its  database,  while  the  IDeX  calculation  used  a  closer 
estimate  of  the  composition  (C70H110O1).  Several  IDeX 
calculations  were  performed  with  this  binder  replaced  by 
typical  plasticisers  —  in  all  cases,  the  results  were  similar, 
with  nothing  approaching  the  BKW  results. 
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Table  2  summarises  die  predicted  equilibnum  detonation 
products  There  are  slight  differences  due  to  the  codes 
allowing  different  product  species.  However,  the  major  effect 
is  believed  to  be  due  to  the  differences  in  the  compress¬ 
ibilities  between  the  Cowan  and  the  Mumaghan  EoS 
descriptions  of  graphite.  This  would  affect  the  CfgnphneVCO 
balance,  causing  a  subsequent  shift  in  the  Hy/HyO  balance. 


TABLE  2.  PREDICTED  DETONATION  PRODUCTS 


Species 

Mole  numbers  (1  kg  explosive)  | 

BKW 

IDeX 

CHa 

n.a.* 

0.407 

CO 

0.443 

5.30 

C02 

00158 

0.232 

H: 

8.66 

13.7 

H 

0  0510 

n.a.* 

nh3 

n.a.* 

1.69  x  10-’ 

HyO 

6.75 

0.498 

HQ 

2.90 

3.66 

Clj 

0.382 

Nil 

n2 

4.53 

4.53 

NO 

4.09  x  10-J 

Nil 

C  (diamond) 

n.a.# 

Nil 

C  (graphite) 

9.64 

5.46 

AlyOy 

4.63 

4.63 

*  Not  available  in  product  species  database. 


The  IDeX  predictions  have  been  adopted,  for  the 
pragmatic  reason  that  the  subsequent  CPeX  analysis 
(described  below)  based  on  the  BKW  predictions  could  not 
be  made  to  fit  the  experimental  detonation  velocity  data. 

REACTIVE  EQUATION  OF  STATE  MODEL 

The  unreacted  or  explosive  phase  is  described  by  a  Mie- 
Grtineisen  EoS  in  the  form 


The  reacted  or  product  phase,  denote  1  by  the  subsenpt  p, 
is  described  by  a  polytropic  EoS  with  a  density-dependent 
index,  namely 


e 


Ppvp 


(7) 


where  Yp  =  ro*  fi/v,  +  y:/vp 


(8) 


and  where  qp  is  the  heat  of  reaction  The  constants  in 
equation  (8)  are  determined  by  requiring  that  it  return  the 
correct  values  for  (t?  In  p/<?  Ln  v)^  at  the  ideal  CJ  stale  and  at 
infinite  expansion. 


The  equation  of  state  for  the  reacting  mixture  is  then 
completely  specified  by  invoking  the  simple  mixture  rules 

P  *  Pi  a  Pp  • v  *  vi  3  vp  ■  t  =  (1  -  A)e.,  +  A  ep  (9) 

where  A  is  the  extent  of  reaction,  varying  from  0  for  the 
unreacted  explosive  to  1  for  the  detonation  products.  The 
associated  reaction  rate  was  developed  for  composite  porous 
explosives  by  Kirby  and  Leiper,9  namely 


A.(i-A)fa«fS.  +  £fi.  +  £fLl 

1  T-  V  I 

here 

[P-Pcf  for  P>  4pCf/3 


GO) 


(U) 


The  subscripts  are:  h  for  hotspot,  i  for  intermediate,  and / 
for  final  stages  of  the  reaction.  There  are  four  adjustable 
parameters  —  three  characteristic  reaction  times  r,  and  the 
critical  pressure  pcr  that  inhibits  the  onset  of  the  hotspot 
reaction.  The  a  factors  in  equation  (10)  describe  the  assumed 
geometry  of  the  bum  front,  controlling  the  switching  on  and 
off  of  the  hotspot,  intermediate  and  final  reaction  rate  terras. 
They  are  functions  of  A,  and  are  Gaussian  in  shape,  namely 


e.  =  e, 


+  ^(P*-Pr) 


(1) 


where  p  is  pressure,  v  is  specific  volume,  e  is  specific  internal 
energy,  and  1~  is  the  Griineisen  coefficient.  The  subscripts 
are:  x  for  the  unreacted  explosive,  0  for  the  initial  state,  and  r 
for  the  reference  (principal)  isentrope.  The  latter  is 
represented  by  the  Birch-Mumaghan  finite  strain  equation  in 
the  form  described  by  Jeanloz8  as 


and 


where 


Pr  =3Koe>(29  +  l)5/-(l  +  <J,ip] 
«r=!*OlVo  9>2(l+(3<3i/3)<»] 
P  =  ^[(vio/wx)2/5-l] 


(2) 

(3) 

(4) 


The  required  coefficients  in  equations  (2)  to  (4)  are 
determined  from  the  Hugoniot 

U  =  cq+su  (5) 

where  U  is  shock  velocity  and  u  particle  velocity,  as 

*0=Px0c0  *nd  ai=(V2K4r-5)  (6) 


exp1 

exp 


-[(a-c,)/^] 

-[(A-c.yw.f 


for  0  S  A  S  C, 
for  Cj  <  A  S 1 


(12) 


0  for  0  S  A  S  C, 

l-a*  for  C,  <  A  S  C/  (13) 

exp|— [( A  -  C,  )/wh  j2  j  -  a*  forC/<A<l 


0 

'"{l-flA-O, 


for  OS  A  SC/ 
for  Cf  <  A  S 1 


(14) 


where  the  Gaussian  parameters  are  defined  in  terms  of  the 
mass  fractions.  O,  of  the  three  stages  as: 


Centroids 


Half-widths 


c,=<t>; 

C/=(o*+o1): 
W*=2C,/(l  +  C,) 
*,=*/(!-«,) 


(15) 


(16) 


'*/ 


k 


•  • 


•  •  • 


•  • 


•  • 


•  • 
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STEADY  NONIDEAL  DETONATION  -  CPe.X  ANALYSIS 

The  above  reactive  equation  of  stale  was  embedded  into 
the  analy  tical  non-ideaj  detonation  model  CPeX  (standing  for 
Commercial  Eerfomunce  of  explosives).  This  model  was 
developed  to  describe  the  detonics  of  non  ideal  explosives  by 
Kirby  and  Leiper’  who  extended  the  small  divergent 
detonation  theory  of  Wood  and  Kirkwood.10  It  describes  the 
flow  along  the  central  streamtube  between  the  detonation 
front  and  the  CJ  plane  for  unconflned  cylindrical  geometry  by 
using  the  shooting  method  described  by  Braithwaite  et  al.n  to 
solve  the  system  of  partial  differential  equations  including  the 
Euler  equations  of  motion,  the  equation  of  state,  the  reaction 
rate,  and  an  empirical  relationship  between  the  wave  front 
curvature  Rs  and  the  charge  diameter  d.  In  previous  versions 
of  CPeX.  this  relationship  had  the  form 

R,  =(d-  Bxqj  )/ A  ^7) 

where  xcj  is  the  CJ  zone  length  at  the  diameter  d,  and  where 
A  and  B  are  empirical  constants.  However,  in  practice,  xqj  is 
not  known  beforehand,  so  that  an  iterative  procedure  must  be 
employed  to  determine  the  necessary  model  constants.  This 
work  employs  a  simpler  form  suggested  by  Leiper,*2  namely 

±  =  a  +  (18) 

Rs  d 

where  a  and  P  are  empirical  constants  and  dcr  is  the  critical 
diameter  below  which  detonation  will  not  propagate. 

The  adjustable  parameters  in  the  reaction  rate  law  were 
varied  until  the  predicted  variation  of  detonation  velocity 
with  unconflned  charge  diameter  matched  the  experimental 
measurements  of  Forbes  et  al.1  as  shown  in  Figure  1.  This 
figure  also  includes  CPeX  fits  to  an  ideal  and  a  non-ideal 
explosive  for  comparative  purposes.  The  former  is  the  fine 
RDX  Composition  B  (COMP  B)  data  of  Malin.13  The  latter 
is  a  proprietary  ICI  blend  of  porous  ammonium  nitrate  prill 
plus  fuel  oil  (ANFO)  with  an  AN/water-in-oil  emulsion. 


Cntical  tfamatw  (d^  /  Char?#  jiimtiir 


FIGURE  1.  DETONATION  VELOCITY  FOR 
UNCONFINED  CYLINDRICAL  CH.ARGES 

The  ideal  detonation  velocity  assumed  in  the  CPeX  fit  for 
PBXW-1 15  is  about  0.5  mm.us'1  faster  than  that  determined 
by  a  linear  extrapolation  of  Forbes'  data.  It  is  believed  that 
the  experimental  data  have  been  obtained  over  too  small  a 
range  of  diameters,  being  limited  to  within  a  factor  of  two  of 
the  critical  diameter.  As  an  illustrative  example,  the  "ideal" 
detonation  velocity  inferred  for  HANFO  from  measurements 
made  on  charges  up  to  twice  the  critical  diameter  is  about 
1.5  mm.ps'1  lower  than  the  value  predicted  from  ideal 


thermodynamic  calculations,  and  inferred  from  measurements 
over  a  wider  range  of  charge  diameters  (namely  up  to  six 
times  the  cntical  diameter). 

The  parameter  values  for  PBXW-115  are  summarised  in 
Table  3.  applying  specifically  to  the  NSVVC  composition.1 : 
For  reasons  not  yet  identified,  the  Australian  variant4  has  a 
cntical  diameter  roughly  twice  that  of  the  NSWC  material 

Since  a  Hugomot  has  not  been  published  for  PBXW-1 15. 
one  was  denved  from  those  of  its  components  using  the 
mixture  rule  desenbed  by  Afanasenkov  et  al.14  —  this 
Hugoniot  is  similar  to  that  measured  experimentally  for  the 
propellant  SP1S-44  with  composition  20/49/21/10% 
HMX  /  AP  /  Al  /  binder,  for  which  U  =  2.774  +  1  855u  15 


TABLE  3.  CPeX  MODEL  CONSTANTS  FOR  PBXW-1 15 


Parameter 

Symbol 

Value 

Initial  density 

PC 

1.79  g  cm  5 

Hugoniot  intercept 

CO 

2.80  mm  |as' 1 

Hugoniot  slope 

s 

1.83 

Griineisen  coefficient 

Txo 

2.60 

Curvature  parameter 

a 

0.165 

Curvature  parameter 

p 

0.692 

Polytropic  coefficient 

X) 

1.343 

Polytropic  coefficient 

n 

0  2045  cm3  g-* 

Polytropic  coefficient 

n 

0.005112  cm6  g  : 

Hotspot  mass  fraction 

0.15  i 

Intermediate  mass  fraction 

0.60 

Final  mass  fraction 

<t>f 

0.25 

Critical  hotspot  pressure 

Per 

5.5  GPa 

Hotspot  time  constant 

T ft 

ll.8jts.CPa-* 

Intermediate  time  constant 

ti 

100  Ops.  GPa  * 

Final  time  constant 

V 

66  0  ps.  GPa-1 

Initially,  the  CPeX  fit  for  PBXW-115  was  based  on  the 
plausible  assumption  that  the  hotspot  or  initial  reaction 
consumed  the  RDX.  that  the  intermediate  phase  of  the 
reaction  involved  the  AP  plus  binder,  and  that  the  Al  was  the 
last  to  react.  The  appropriate  mass  fractions  would  then  be 
<&h  =  20%,  =  55%  and  <ly=  25%.  However,  a  superior  fit 

to  the  velocity  data  was  achieved  by  modifying  the  hotspot 
mass  fraction  to  be  <t>&  =  15%  (with  <t>,  =  60%.)  This  may 
indicate  that  the  Gaussian  shape  factors  developed  for  highly 
porous  explosives  and  described  in  equations  12  to  16  are  not 
as  appropriate  for  non-porous  explosives.  Alternatively,  it 
may  be  argued  that  the  use  of  a  single  parameter,  A,  is  not 
sufficient  to  simultaneously  account  for  the  release  of  energy 
via  the  chemical  reactions,  the  switch  from  unreacted  to 
reacted  equations  of  state,  and  the  geometry  of  the  bum  front 

It  should  be  noted  that  the  adopted  values  of  the  constants 
describing  the  late  stages  of  the  reaction  (ie.  dy,  t,  and  trt 
are  not  unique  —  other  assumed  values  could  give  equally 
acceptable  fits  to  the  detonation  velocity  data.  This  is 
explained  by  examining  the  reaction  rate  in  greater  detail. 
Figure  2  shows  how  the  predicted  C7  extent  of  reaction  varies 
with  charge  diameter.  Detonation  in  the  ideal  explosive 
COMP  B  fails  when  its  CJ  extent  of  reaction  falls  below 
96%,  whereas  both  the  non-ideal  compositions  continue  to 
propagate  down  to  about  15%  reaction.  The  maximum  CJ 
extent  of  reaction  achieved  over  the  range  of  diameters 
sampled  for  PBXW-115  is  about  35%.  Consequently,  that 
part  of  the  reaction  rate  for  A  >0.35  has  not  been  fitted  to 
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FIGURE  2.  CPeX  CJ  EXTENT  OF  REACTION 


experimental  data,  and  so  must  be  considered  conjectural. 

Figure  3  examines  this  in  greater  detail,  by  presenting  the 
full  reaction  rate  surface  X(p.  A)  for  each  explosive.  The 
(pcj-Xcj)  locus  (traced  out  by  varying  charge  diameter)  for 
unconfined  detonation  is  superimposed  on  this  surface, 
together  with  the  two  (p.A)  histories  followed  by  Lagrangian 


particles  on  the  cental  streamtube  between  the  shock  front 
and  the  CJ  plane  for  -neon fined  detonations  in  the  smallest 
and  largest  diameters  cturaciensed  experimentally  It  con  be 
seen  that  the  detonation  velocity  measurements  performed  to 
date  have  sampled  only  a  very  restneted  subset  of  the 
reaction  rate  surface  for  PBXW  1 15,  so  that  any  predictions 
of  late  time  reactions  involving  the  AP  and  Al  must  be  of 
dubious  accuracy  In  order  to  extract  accurate  values  for  the 
time  constants  t,  and  yfrom  a  CPeX  analysis,  it  would  be 
necessary  to  measure  tne  detonation  velocity  of  unconfined 
charges  exceeding  250  mm  diameter  and  1500  mm  length. 

The  measurements  for  COMP  B  have  sampled  across  the 
full  range  of  extents  of  reaction,  but  only  for  a  restricted  part 
of  the  pressure  domain.  Consequently,  descriptions  of  low 
pressure  shock  initiation  based  only  on  the  CPeX  fit  would  be 
suspect  unless  augmented  by  specific  initiation  experiments. 
However,  the  reaction  rate  surface  for  HANFO  is  very  well 
covered  by  the  experimental  detonation  velocity  data, 
permitting  accuate  predictions  of  its  detonics. 

CJ PRESSURE 
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FIGURE  3.  CONTOURS  OF  REACTION  RATE  A(p.A) 


Figure  4  presents  the  CPeX  predictions  for  the  pressures 
at  the  shock  front,  and  at  the  CJ  plane.  It  can  be  seen  that  they 
are  predicted  to  vary  by  a  factor  of  about  three  over  the  full 
range  of  diameters  that  will  support  detonation. 


FIGURE  4.  CPeX  SHOCK  AND  CJ  PRESSURES 

HYDROCODE  MODELLING  OF  PBXW-l  15 

The  equation  of  state  and  reaction  rate  law  described 
above  for  the  CPeX  model  were  embedded  in  the  explicit 
finite  element  hydrocode  DYNA2D3  to  permit  the  simulation 
of  time-dependent  reactive  flow  in  non-ideal  explosives.  The 
constants  derived  for  PBXW-l  15  by  the  CPeX  analysis,  and 
displayed  in  Table  3.  were  used  in  the  DYNA2D  simulations 
without  any  further  adjustment. 

STEADY-STATE  AXISYMMETRIC  DETONATION 

Figure  5  summarises  the  DYNA2D  predictions  of 
detonation  velocity  in  axisymmetric  geometry.  AU  charges 
were  3C0  mm  in  length  to  ensure  ih3t  steady  state  conditions 
were  reached.  There  is  excellent  agreement  between  the 
numerical  simulations  and  the  experimental  data1  for  both 
unconfined  and  confined  detonations  in  2.5  mm  thick  brass. 

The  experimental1  pass  /fail  diameters  for  PBXW-115 
were  38.7  /  35.1  mm  unconfined.  and  22.2  / 19  mm  confined 
(in  2.5  mm  brass).  The  simulations  were  in  excellent 
agreement,  predicting  38  /  37  mm  uncon  fined  and  22/19  mm 
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FIGURE  5.  PBXW-1 15  DETONATION  VELOCITY 

confined.  Figure  6  shows  how  the  predicted  shock  velocity  in 
a  37  mm  rate  stick  decayed  with  distance  from  the  booster. 
The  wave  propagated  for  over  four  charge  diameters  before 
dying.  Dick14  has  observed  similar  behaviour  for  a  HMX- 
based  composite  propellant.  This  suggests  that  the  153  mm 
length  charges  used  at  NSWC  were  only  just  long  enough  to 
ensure  steady  state  behaviour  in  diameters  close  to  critical. 


FIGURE  6.  FAILURE  IN  3 7 MM  DIAMETER 
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measurements1  and  the  DYNA2D  simulations  of  shock  front 
radius  of  curvature  for  unconfined  charges,  plotted  to 
illustrate  the  linear  behaviour  expressed  in  equation  (18).  The 
scaner  reflects  the  difficulty  in  measuring  the  curvature  both 
experimentally  and  in  the  simulations  This  data  was  used  in 
the  CPeX  fit.  but  not  explicitly  by  the  DYNA2D  simulations. 
The  good  agreement  thus  helps  to  confirm  the  flow  pattern 
assumed  in  the  CPeX  theory. 

Wave  curvature  experiments  are  often  performed  in  order 
to  estimate  the  CJ  zone  lengths  via  application  of  the  Wood- 
Kirkwood  theory.10  The  zone  lengths  calculated  by  Forbes  et 
al.1  are  compared  with  the  CPeX  predictions  in  Figure  8.  The 
agreement  is  poor,  with  the  CPeX  CJ  zone  lengths  being  up 
to  a  factor  of  two  higher  than  Forbes'  calculations  based  on 

?-:d.=c!zl 


where  D«  is  the  ideal  value  of  the  detonation  velocity,  D  the 
measured  value,  and  Rs  the  wave  curvature  for  the  charge 
diameter  of  interest.  C  is  a  constant  dependent  upon  the 
equation  of  state.  Wood  and  Kirkwood  derived  equation  (19) 
as  a  specific  solution  to  their  more  general  theory  by  making 
a  number  of  assumptions.  In  particular,  they  assumed  that  the 
reaction  rate  is  Arrhenius  in  form,  so  that  the  reaction  rate 
immediately  behind  the  shock  is  low.  and  then  rises  to  a  sharp 
maximum  very  close  to  complete  reaction.  This  leads  firstly 
to  CJ  parameters  (p .  v  and  A),  that  are  only  slightly  perturbed 
from  their  ideal  values,  and  secondly  to  a  pressure  profile  that 
is  essentially  fiat-topped  at  the  von  Neumann  spike  value. 
These  conditions  are  not  too  dissimilar  from  the  behaviour 
summarised  in  Figures  2  and  3  for  the  ideal  explosive  COMP 
B.  However,  these  Figures  show  that  these  assumptions  are 
clearly  inappropriate  for  both  non-ideal  explosives.  PBXW- 
115  and  HANFO.  Hence,  the  simplistic  form  of  the  Wood- 
Kirk  wood  theory,  represented  by  equation  (19),  should  not  be 
applied  to  non-ideal  explosives. 


This  high  level  of  agreement  for  unconfined  detonation 
confirms  not  only  the  model  constants,  but  also  serves  to 
verify  the  assumptions  made  in  the  CPeX  model  to  relate 
flow  divergence  to  shock  front  curvature.  Such  assumptions 
are  not  required  for  the  DYNA2D  simulations.  Furthermore, 
the  similar  excellent  agreement  for  confined  detonation 
provides  additional  confirmation  of  the  model,  having 
sampled  a  different  combination  of  lateral  expansion  and 
extent  of  reaction. 

WAVE  FRONT  CURVATURE  AND  CJ  ZONE  LENGTHS 
Figure  7  presents  the  comparison  between  the  NSWC 


FIGURE  7  SHOCK  FRONT  CURVATURE 


FIGURE  8.  CPeX  CJ  ZONE  LENGTHS 


It  is  apparent  from  Figure  8  that  the  CJ  zones  for  the  two 
non-ideal  explosives  extend  over  much  larger  regions  of 
space  than  is  the  case  for  the  ideal  one.  In  particular,  it  is 
predicted  that  experiments  intended  to  characterise  an  ideal 
steady-state  planar  detonation  in  PBXW-1 15  must  be 
performed  on  a  massive  scale  —  the  ideal  CJ  zone  predicted 
by  CPeX  exceeds  60  mm  in  length  and  15  ps  in  duration. 

CORNER  TURNING 

Forbes  et  al.'-2  have  characterised  the  comer-turning 
ability  of  PBXW-1  IS.  measuring  the  breakout  times  through 
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the  curved  surfaces  of  bowl  shaped  acceptor  charges  The 
bowls  were  initiated  through  their  flat  rear  surfaces  by 
boosters  of  51  mm  diameter  and  153  mm  in  length.  The 
curved  surface  of  each  bowl  had  a  radius  of  curvature  of 
51  mm  (measured  from  the  edge  of  the  booster)  Four 
different  booster  configurations  were  used  —  (l)  bare 
PB.XW-l  15.  <1 )  bare  COMP  B,  (3)  bare  COMP  B  with  a  steel 
plate  covering  the  flat  surface  of  the  bowl,  and  (4)  PBXW- 
1 15  confined  in  a  brass  tube  with  16  5  mm  thick  walls. 

In  order  to  control  the  distortion  of  the  finite  element 
mesh  at  the  booster/bowl  junction,  it  was  necessary  to 
introduce  a  6  mm  radius  of  curvature  there  in  the  numerical 
simulations  with  DYNA2D  —  the  effect  of  this  change  on  the 
predicted  breakout  times  is  unknown,  but  is  believed  to  be 
small  Figure  9a  presents  the  predicted  pressure  contours  9  |is 
after  the  shock  from  the  booster  has  entered  the  bowl,  while 
Figure  9b  shown  the  extent  of  reaction  after  the  shock  has 
completely  enveloped  the  bowl  (with  the  material  motion 
removed  in  order  to  facilitate  comparison). 

The  comparison  with  experiment  is  shown  in  Figure  10, 
and  is  seen  to  be  excellent  for  the  bare  boosters,  cases  (1), 
(2),  and  (3).  This  agreement  indicates  that  the  reaction  rate 
surface  for  PBXW-115  (shown  in  Figure  3a)  is  accurate  in 
regions  beyond  that  directly  calibrated  by  the  detonation 
velocity  measurements.  In  particular,  the  breakout  through 
the  outer  edges  of  the  bowl  for  case  (1),  the  unconfined 
PBXW-1 15  booster,  is  controlled  by  the  reaction  rate  surface 
in  the  lower  pressure  regime  below  about  5  CPa.  The  break¬ 
outs  for  cases  (2)  and  (3)  with  the  COMP  B  boosters  are 
dominated  by  the  high  pressures,  exceeding  20  GPa, 
developed  along  the  axis  of  the  bowls. 

The  lack  of  agreement  for  case  (4).  the  confined  PBXW- 
1 15  booster,  was  unexpected.  Both  the  experimental  detona¬ 
tion  velocity  data  and  the  associated  simulations  summarised 
in  Figure  5  prove  that  even  thin  (2.5  mm)  brass  confinement 
makes  a  significant  difference  to  detonation  in  PBXW-115, 
decreasing  the  critical  diameter  by  a  factor  of  roughly  two. 
The  confinement  in  the  comer  turning  experiment  was  a  brass 
tube  with  thick  (16.5  mm)  walls,  and  hence  should  increase 
both  the  peak  value  and  the  duration  of  the  pressure  delivered 


FIGURE  10  CORNER  TURNING  BREAKOUT  TIMES 


by  the  booster  to  the  acceptor  charge  Furthermore,  the  shock 
wave  in  the  brass  confinement  enters  the  bowl  in  concert  with 
that  from  the  booster,  increasing  the  area  over  which  pressure 
is  applied  to  the  bowl.  These  factors  should  all  combine  to 
appreciably  shorten  the  breakout  times  through  the  outer 
edges  of  the  bowl.  No  plausible  explanation  for  the  discrep¬ 
ancy  between  experiment  and  simulation  could  be  found. 

WEDGE  TEST  AND  POP-PLOTS 

The  wedge  test  is  often  used  to  characterise  the  shock 
initiation  behaviour  of  explosives.  A  plane-wave  generator 
sends  a  relatively  flat-topped  planar  shock  into  a  wedge  of 
acceptor  explosive  —  observation  of  the  emergence  of  the 
shock  through  the  angled  face  of  the  wedge  gives  a 
continuous  record  of  run  up  to  detonation.  The  results  are 
displayed  on  a  Pop-plot,  where  the  measured  run  distance  to 
detonation  is  plotted  against  the  initial  shock  pressure. 

Although  PBXW-115  has  not  been  investigated 
experimentally  by  this  technique,  several  simulations  were 
performed  to  predict  its  likely  behaviour.  Figure  II 
summarises  the  predicted  shock  velocity-time  trajectories 
over  a  50  mm  run  distance  induced  by  various  input  shock 


(i)  sv«  oaxw-ns  (2)  airvccup  a  <oi  comp  3'  SwipWit  («icoor«t<jpsxw-it3 


(a)  Pressure  after  shock  enters  bowt 


(b)  Extent  of  reacton  after  shock  has  oompieiety  enveloped  bowl. 
(The  distortion  has  boon  romovod  to  faciitate  comparison ) 


FIGURE  9  EFFECT  OF  BOOSTER  TYPE  ON  CORNER  TURNING  LN  BOWL  OF  PBXW-115. 
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FIGURE  11.  PREDICTED  SHOCK  ACCELERATION 
OVER  50  MM  IN  WEDGE  TEST 

pressures.  In  the  simulations,  these  pressures  were  created  by 
the  impact  of  thick  aluminium  flyer  plates. 

These  trajectories  are  notable  for  the  absence  of  the  sharp 
transition  between  low  velocity  shock  and  high  velocity 
detonation  that  is  characteristically  observed  with  ideal 
explosives.  They  are  qualitatively  similar  to  those  published 
for  other  non-ideal  explosives.  Dick16  suggested  that  the 
gradual  increase  of  velocity  with  time  shown  by  the  traces  for 
a  HMX  /  AP  /  A1  /  binder  propellant  (not  unlike  PBX  W- 1 15) 
was  evidence  of  overdriven  detonation.  In  spite  of  all 
trajectories  showing  gradual  acceleration,  workers  at 
RCEM17'1*  derived  a  Pop-plot  for  an  emulsion  explosive, 
concluding  that  it  was  much  less  shock  sensitive  than  TNT. 

The  present  work  suggests  that  these  interpretations  of 
wedge  test  data  for  non-ideal  explosives  are  incorrect.  Figure 
12  extends  the  trajectories  of  Figure  1 1  out  to  a  distance  of 
500  mm.  Two  features  are  apparent.  Firstly,  the  trajectories 
for  the  2.2  and  3.6  GPa  initial  shocks  now  show  clear 
transitions  between  low  and  high  velocity  regimes.  Secondly, 
none  of  the  trajectories  have  attained  a  steady-state  velocity 
even  after  500  mm  of  run.  though  they  are  all  asymptoting 
towards  the  expected  ideal  detonation  velocity  of 
6  67  mm. ns1. 


O&anc*  (mm) 

FIGURE  12.  PREDICTED  SHOCK  ACCELERATION 
OVER  500  MM  IN  WEDGE  TEST 

This  behaviour  can  be  explained  by  referring  to  the 
reaction  rate  surfaces  in  Figure  3.  For  PBXW-11S.  the  A 
surface  can  be  divided  into  two  regimes. 

Firstly,  for  initial  shock  pressures  above  5  GPa.  the  peak 
reaction  rate  occurs  at  minimal  extents  of  reaction,  causing 
rapid  consumption  of  the  RDX  component.  Any  rapid 
acceleration  of  the  shock  velocity  will  occur  during  this  short 
phase,  quickly  giving  rise  to  velocities  typical  of  detonation 


in  small  diameters,  re  about  5  0  »  5  5  mm  ns  1  However, 
once  the  extent  of  reaction  exceeds  20 %  or  so.  the  reaction 
rate  decreases  sharply,  and  so  approach  to  the  ideal 
detonation  state  will  be  gradual. 

Secondly,  for  initial  shock  pressures  below  5  GPa.  the 
initial  reaction  rate  is  tow,  so  that  the  shock  propagates  with 
minimal  acceleration.  However,  as  the  reaction  proceeds  and 
the  pressure  builds  up.  the  condition  is  reached  where  the 
reaction  rate  increases  sharply  The  shock  then  makes  a  sharp 
transition  to  a  higher  velocity,  though,  as  above,  it  is 
supported  essentially  by  consumption  of  the  RDX  Further 
acceleration  is  consequently  slow. 

In  contrast,  ideal  explosives  at  low  porosity  will  have 
reaction  rate  surfaces  similar  to  that  shown  in  Figure  3c  for 
COMP  B.  Here,  the  initial  reaction  rate  at  low  to  medium 
shock  pressures  is  low.  so  that  initial  shock  velocities  are  low 
with  minimal  acceleration.  However,  once  several  percent 
reaction  has  occurred,  the  increase  of  pressure  causes  the 
rapid  acceleration  of  the  reaction  rate,  resulting  in  an  abrupt 
transition.  Furthermore,  the  reaction  rate  remains  high,  so  that 
the  reaction  soon  goes  to  completion,  and  the  ideal  detonation 
velocity  is  quickly  attained.  Hence,  there  is  always  a  strong 
distinction  between  detonation  and  non-detonation. 

This  analysis  suggests  that  reduction  of  wedge  test  traces, 
and  interpretation  of  their  resultant  Pop-plots,  requires 
substantial  modification  for  non-ideal  explosives.  Pop-plots 
are  traditionally  used  for  two  distinct  purposes.  Firstly,  they 
are  used  to  derive  kinetics  for  detonation  models,  as 
exemplified  by  the  Forest-Fire  model1’.  However,  this  paper 
has  shown  how  this  can  be  accomplished  instead  by 
detonation  velocity  measurements  on  non-ideal  explosives. 
Secondly,  they  are  used  to  rank  the  shock  sensitivity  of 
different  explosives  for  hazard  studies.  For  this  purpose,  it  is 
sufficient  to  know  the  conditions  under  which  incident  shocks 
lead  to  self-propagating  reactive  waves  releasing  substantial 
amounts  of  energy.  Hence,  it  is  necessary  to  replace  the 
identification  of  transition  to  detonation  by  the  identification 
of  transition  between  low  and  high  velocity  shocks. 

Once  this  is  done,  Figures  11  and  12  show  that  PBXW- 
115  exhibits  such  transitions  for  incident  shock  pressures 
below  5.0  GPa.  The  resulting  Pop-plot  is  shown  in  Figure  13. 
where  it  can  be  seen  that  PBXW-1I5  exhibits  a  greater 
hazard  to  low  strength  shocks  than  either  TNT  or  Baratol. 
both  of  which  have  much  smaller  critical  diameters. 


FIGURE  1 3.  PREDICTED  POP-PLOT  FOR  PBXW-1 15 


The  wedge  test  behaviour  of  other  non-ideal  explosives 
can  now  be  explained. 
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In  the  propellant  srndied  by  Dick14,  HMX  would  play  the 
role  of  the  RDX  m  PBXW115  —  Figure  1  of  Dick's  paper 
does  indeed  show  a  strong  similarity  to  Figure  11  of  this 
work,  though  shifted  to  higher  initial  shock  pressures  in  line 
with  the  lower  inherent  shock  sensitivity  of  HMX  when 
compared  with  RDX 

The  emulsion  studied  at  RCEM17 11  is  expected  to  have  a 
reaction  rate  surface  not  unlike  that  of  the  HANFO  shown  in 
Figure  3b.  This  exhibits  only  the  first  regime  discussed 
above,  where  the  reaction  rate  is  always  greatest  at  or 
immediately  behind  the  shock  front,  so  that  acceleration  is 
always  smooth,  and  transitions  would  never  occur.  Hence, 
both  the  published  Pop-plot  for  this  emulsion,  and  the  claim 
that  it  is  less  sensitive  than  TNT,  are  suspect. 

IDEAL  DETONATION  VELOCITY  EXPERIMENT 

Held4  has  recently  described  a  technique  to  measure  the 
ideal  detonation  velocity  using  only  small  samples,  and  has 
applied  it  to  the  German  variant  of  PBXW-115.  A  donor 
charge  with  high  detonation  velocity  drives  a  bow  wave  into 
a  parallel  adjacent  acceptor  test  charge  with  lower  detonation 
velocity,  which  is  inferred  from  the  measurement  of  the 
breakout  through  the  face  opposite  the  initiation  point.  Held 
assumed  that  if  a  steady  detonation  was  established  in  the 


(b)  Extent  of  Reaction 

FIGURE  14.  DYNA2D  PROFILES  FOR  IDEAL 
DETONATION  VELOCITY  METHOD 


acceptor,  then  it  must  represent  an  ideal  detonation. 

Figure  14  shows  the  predicted  detonation  wave  shapes 
immediately  pnor  to  breakout  for  the  case  of  an  Ociol  donor 
dnving  a  bow  wave  into  PBXW- 1 15.  Although  the  bow  wave 
in  the  PBXW-115  is  shown  as  being  curved,  DYNA2D 
predicts  that  it  does  become  straight  prior  to  full  breakout 
(This  simulation  had  to  be  performed  using  slab  charges  in 
planar  geometry,  in  place  of  the  actual  half  cylinders.  The 
PBXW-115  used  in  the  experiments  had  the  composition 
40/24/  24/1 2%  AP  /  RDX  /  A1  /  binder.  This  was  treated 
by  increasing  to  1ST  and  decreasing  <t>,  to  57%  ) 

The  charges  used  by  Held  were  60  mm  in  length,  which  is 
roughly  the  same  order  as  the  ideal  CJ  zone  length  predicted 
by  CPeX  (shown  in  Figure  8)  Consequently,  DYNA2D  is 
predicting  that  a  steady  non-ideal  detonation  has  been  created 
in  the  PBXW- 1 15.  with  the  extent  of  reaction  having  reached 
a  maximum  of  less  than  by  50%  by  the  time  that  the 
detonation  from  begins  to  break  out  through  the  end  face. 

The  predicted  breakout  times  for  the  non-ideal  detonation 
are  compared  with  the  experimental  results  in  Figure  15.  with 
the  agreement  being  excellent,  confirming  the  accuracy  of  the 
DYNA2D  simulation.  The  conclusion  is  that  the  small 
sample  technique  does  not  measure  ideal  detonation  velocity 
as  intended.  Much  longer  and  larger  diameter  charges  would 
be  required  to  accomplish  this  objective. 


Rad'us  (mm) 

FIGURE  15.  BREAKOUT  TIMES  FOR  IDEAL 
DETONATION  VELOCITY  METHOD 


CONCLUDING  REMARKS 

The  application  of  the  CPeX  small-divergent  detonation 
theory,  and  the  DYNA2D  simulations  performed  for  this 
paper,  have  both  been  built  around  the  assumption  that  the 
ideal  detonation  velocity  of  PBXW- 1 15  is  6.67  mm.ps'1  as 
predicted  by  the  chemical  equilibrium  code  IDeX.7  At  the 
time  of  writing,  there  is  no  direct  experimental  evidence  to 
support  this  assumption.  However,  the  ensuing  excellent 
agreement  between  the  predictions  and  almost  all  the 
available  experimental  data  provides  strong  indirect  evidence 
that  this  assumption  is  correct,  and  hence  that  the  deionics  of 
PBXW-1 15  are  strongly  non-ideal. 

The  only  piece  of  experimental  data  which  could  not  be 
reproduced  by  the  DYNA2D  simulations  was  the  breakout 
time  in  one  configuration  of  the  comer  turning  experiments. 
Strangely,  the  simulations  did  successfully  reproduce  the 
breakout  limes  of  the  two  configurations  that  were  found 
experimentally  to  bracket  the  offending  configuration. 

The  analysis  presented  in  this  paper  illustrates  that  the 
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reaction  rate  surface  of  the  non  ideal  explosive  PBXW-1 15  is 
fundamentally  different  from  that  of  the  ideal  explosive. 
Composition  B  It  is  strongly  believed  that  this  is  a  specific 
example  of  a  more  general  truth 

Many  of  the  experimental  tests  in  use  today  were 
developed  to  characterise  the  detorucs  of  ideal  explosives, 
and  at  a  fundamental  level,  require  the  reaction  rate  surface  to 
have  particular  features  When  applying  such  tests  to  non¬ 
ideal  compositions,  it  is  essential  to  understand  the  basis  of 
what  the  test  is  actually  measuring  in  order  to  know  if  it  is 
still  applicable.  In  particular,  if  knowledge  of  the  full  range  of 
detonation  behaviour  is  required  for  an  explosive  like 
P8XW-I15,  there  is  no  way  to  avoid  at  least  some 
experiments  performed  on  massive  charges. 

Many  of  the  theories  and  the  data  reduction  techniques 
that  are  applied  to  detonation  experiments  assume  that  non¬ 
idea!  behaviour  can  be  treated  as  a  small  perturbation  in  (p.  v. 
A.)  space  from  the  ideal  CJ  conditions  However,  a  steady- 
state  detonation  that  is  supported  by  only  about  I5?c  reaction 
cannot  be  considered  to  represent  a  small  perturbation  from 
ideal.  Due  to  the  work  by  Forbes  and  co-workers  at 
NSWC,'-  there  is  now  an  extensive  database  of  experimental 
information  available  on  at  least  one  non-ideal  explosive  for 
use  in  confirming  the  miriad  emerging  theories  of  non-ideal 
detonics. 
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